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Different methods in the study of anomalous scaling of factorial moments in high energy 
e"'"e~collisions are examined in some detail. It is shown that the horizontal and vertical facto- 
rial moments are equivalent only when they are used in combination with the cumulant variables. 
The influence of different reference frames and that of phase space restrictions is also discussed. 
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PACS numbers: 13.85Hd 



I. INTRODUCTION 



Since the observation of exotic multiparticle events in 
Cosmic Ray experiments ||^ and especially the discovery 
of unexpectedly large local fluctuations recorded by the 
JACEE collaboration the investigation of nonlinear 
phenomena (NLP) in high energy collisions has attracted 
much attention [b| . One of the signals of these NLP is the 
fractal property ^ of the multiparticle final states in high 
energy collisions. This property can be characterized by 
the anomalous scaling of the probability moments defined 
as 
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where a certain phase space region A is divided in a 
proper way (isotropically for a self-similar fractal while 
anisotropically for a self-affine fractal into M sub- 
cells, ipra is the probability for a particle to fall into the 
TTith sub-cell, (• • •) denotes the average over an event sam- 
ple. If the Cq(A/)'s have a power law behaviour with the 
diminishing of phase space scale: 



C,(A/) - A/-^' (Af->oo), 



(2) 



then it is said to be anomalous scaling and the system is 
a fractal. 

In real experiments the probability is unobservable 
and the corresponding moments Cq is unaccesible. This 
problem has been solved by Bialas and Peschanski 
who are able to show that the factorial moments 
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are equal to the probability moments Cq provided the 
statistical fluctuations are Poissonian. Thus the scaling 



property of factorial moments, sometimes called intermit- 
tency, becomes a central problem in the study of nonlin- 
ear phenomena in high energy collisions. 

Various methods have been developed in this study. 

Firstly, people noticed that in the definition Eq.(3) of 
factorial moments the average over event sample (vertical 
average) is carried out first and then is the average over 
the M sub-cells (horizontal average) . It was proposed to 
exchange the order of these two averages and define the 
horizontal factorial moments (HFM) ||^ as 
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Accordingly, the Fq defined in Eq.(3) is called vertical 
factorial moments (VFM). 

Note that the equality of factorial moments Fq and 
probability moments Cq has been proved only for the 
VFM. Therefore, in the study of the nonlinear phenom- 
ena — fractal property of multiparticle system, the HFM 
is appropriate only when it is equal to VFM. We will see 
in the following that this equality holds in some cases but 
does not hold in some other cases. 

Secondly, various methods have been proposed to cor- 
rect for the unflatness of the phase-space variable dis- 
tributions. One is to divide the factorial moments by a 
factor Ra 
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Another one is to change the phase space variable x into 
the corresponding cumulant variable Xc before calculat- 
ing the factorial moments. The cumulant variable is de- 
fined as H] 
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FIG. 1: The rotated coordinate system 
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FIG. 2: Average distribution of phase space variables 

where x^in and Xmax are the two boundaries of the x 
region, respectively. 

Another problem arises while carrying on this kind of 
study in e+e^collisions. In these collisions the thrust (or 
sphericity) axis is chosen as the z axis (longitudinal axis) 
to define the phase space variables: rapidity y, transverse 
momentum pt and azimuthal angle ip. Different frames 
could be used to define the azimuthal angle ip. The first 
one is to choose the minor axis of thrust (or sphericity) 
analysis as the x axis, and use it as the starting point 
for counting the azimuthal angle ip. The second one is 
to put the z axis still on the major thrust axis, but turn 
the coordinate system around it and let the new x axis 
lie on the zq-z plan[p^, where xo,yo,ZQ denote the axes 
of the lab system and x, y, z those of the turned system, 
as shown in Fig.l. In the following this frame will be 
referred to as the "rotated frame" . The third method 
is to rotate the frame in each event for a random angle 
around the z axis This is called the random frame. 

All these methods have been used in the literature 
for studying the anomalous scaling of factorial moments, 
making the results hard to be compared. In the present 
paper we will examine these methods in some detail and 
discuss their applicability in physical problems. We will 
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FIG. 3: The distribution of cumulant variables 

take e+e~collisions at the Z° energy y/s ~ 91.2 GeV as 
example and use JETSET7.4 Monte Carlo code to 
generate 500 000 multihadron events as the event sample. 

II. AVERAGE DISTRIBUTIONS OF PHASE 
SPACE VARIABLES 

In Fig. 2 (a), (6) and (d) are shown the average dis- 
tributions of y, Pt and 93, where the rapidity is defined 
as 2/ = 0.5 \n[{E + pz)/{E ~ p^)] with z along the thrust 
axis; the azimuthal angle ip is defined in the plane per- 
pendicular to the thrust axis, calculated with respect to 
the minor axis. 

It can be seen from the figures that all the distributions 
are unflat. Especially, the distribution of pt is exponen- 
tial and is highly concentrated in low pt . A simple vari- 
able transformation p^ — > Inpt psf can make it a litter 
flatter as shown in Fig. 2 (c). 

III. THE CORRECTION FACTOR AND 
CUMULANT VARIABLES 

The unflat average distribution will cause additional 
variation of factorial moments with the diminishing of 
phase space scale and make the scaling property of facto- 
rial moments unequal to that of the probability moments 
even for the VFM. This effect has to be corrected. 

Fialkowski proposed a factor cf. Eq.(5), to correct 
for this effect. This method works good when the dis- 
tribution of the phase space variable is not far from flat, 
e.g. the distribution of rapidity y in a restricted central 
region \y\ < Yc with Fc = 2 as shown in Fig. 2 (a), and 
is not good when the distribution is far from flat. This 
is especially the case for the distribution of pt, cf. Fig. 2 
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ETSET7,4, Correction factor, -2 < y < 7, 0.1 < Pt < 2 
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FIG. 4: VFM (solid circles) and HFM (solid curves) using correction factor method 
(restricted phase space). In (6) the upward triangles are for q — 5, downward ones for q — 4. 




FIG. 5: VFM (solid circles) and HFM (solid curves) using cumulant variables (restricted phase space) 



(b). Therefore, people sometimes transform pt — > Inpt 
first 1 13 and calculate Fq{\npt) instead of Fq{pt), and 
then correct the result by the factor Rq given in Eq.(5). 
Note that the highly asymetric region 0.1 < pt < 2 is 
transformed to the region —2.3 < Inpt < 0.69, which 
is distributed symmetrically around the pick of distribu- 
tion, cf. the two arrows in Fig's. 2 (6) and (c). 

Since a transformation to a flatter distribution is nec- 
essary before calculating factorial moments, it is evident 
that the best way is to transform all the phase space vari- 
ables to a flat distribution first. This could be established 
through the transformation to cumulant variables , cf. 



Eq.(6). The corresponding distributions are shown in 
Fig.3. 



IV. VERTICAL AND HORIZONTAL 
FACTORIAL MOMENTS 

Now let us turn to the comparison of vertical and hor- 
izontal factorial moments (VFM and HFM). 

As noticed in the Introduction, our aim is to explore 
the anomalous scaling of probability moments as shown 
in Eq.(2) but the equivalence of factorial and probability 
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ETSET7,4, Cumulcnt variable, —5 < y < 5, 0,1 < p- < 3 
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FIG. 6: VFM (solid circles) and HFM (solid curves) using correction factor method (nearly full phase space) 



moments has been proved only for the VFM. So, the 
HFM is appropriate only when it is equal to VFM. 

In Figures 4 and 5 (a), (b), (d), (e), (/) are shown the 
1-D and 3-D VFM (sohd circles) and HFM (solid curves) 
for the moment orders g = 2,3,4, 5 calculated using the 
correction factor, Eq.(5), method and the cumulant vari- 
ables Eq.(6), respectively. 

It can be seen from the figures that the VFM and HFM 
are equal only when using together with the cumulant 
variables and are unequal, especially for 1-D Fq{ip) and 
Fq{pt), when using the correction factor method. 

The results from VFM -I- Correction Factor method 
are about the same as that from VFM + Cummulant 
Variables for 1-D Fq{y) and 3-D Fq{y,pt,ip) but are not 
the case for 1-D Fq{ip) and Fq{pt), cf. Fig's. 4 and 5. 



V. THE INFLUENCE OF PHASE SPACE 
REDUCTION 

In hadron-hadron and nucleus-nucleus collisions the 
central rapidity regions with the rapidity y restricted to 
\y\ < for some value of Yc is commonly used. This 
is physically meaningful, because in these collisions the 
final state particles are mainly produced in the central 
region, while the particles in the regions |?;| > Yc are 
mainly come from the fragmentation of incident parti- 
cles (leading particle effect). 

On the contrary, in e+e" collisions, the multihadron fi- 
nal state is produced from a point source — virtual pho- 
ton or Z°, and no leading particle effect is present. The 
rapidity y is usually defined with respect to the thrust or 
sphericity axis. In this case, to carry on the study in a 
"central rapidity region" |y| < 1^ is physically doubtful. 
This is especially evident when the collision energy is so 



high that 3, 4 or even more jets can be produced. In a 
2-jet event the restriction Ij/I < 1^ will cut out the most 
energetic particles from the two opposite jets symmetri- 
cally, but in a 3-jet event the same cut will cut out the 
most energetic particles only from one jet while retain 
almost all the particles in the other two jets. This asym- 
metric cut will results in unexpected phenomena, and the 
physical meaning may be difficult to interpret. Therefore, 
the study of multiparticle dynamics in e"'"e~ collisions is 
better to be carried out in the full phase space. However, 
the central rapidity region is sometimes also used in the 
literature for the study of e+e" collisions Therefore, 
to investigate the influence of rapidity cut is worthwhile. 



In Fig. 6 are shown the results of VFM and HFM in 

nearly the full phase space 5 < y < 5,0 < ip < 

27r,0.1 < pt < 3 GeV, to be compared with the results 

shown in Fig. 5 for a restricted phase space 2 < y < 

2, < < 27r, 0.1 <pt<2 GeV. 



It can be seen from Fig. 6 that the first point in 3-D 
Fq{y,Pt,f) and 1-D Fq{y) and the first 3 points in 1- 
D Fq{ip) do not lie on a scaling curve together with the 
other points. This is due to the momentum conservation 
effect |l4|. In the anomalous-scaling study, these points 
should be omitted. 



The momentum conservation effect will also be reduced 
in a restricted phase space region, which was first pointed 
out in Ref . and has been proposed as a second method 
for eliminating the influence of this effect. This explains 
the reason why the first points in Fig's. 5 lie on the scaling 
curves. 
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VI. THE FLUCTUATIONS IN AZIMUTHAL 
ANGLE 

The fluctuations in azimuthal angle are worthwhile 
special investigation. It is commonly expected that there 
should have cylindrical symmetry around the z axis. If 
that is the case, then the fluctuations should have no cor- 
relation with the X axis chosen for counting the azimuthal 
angle ^p. In Fig's. 4, 5, 6 (6) are shown the Fq{ifi) with 
X axis along the minor of thrust analysis, while in the 
corresponding figures (c) are shown the results after ro- 
tating the X axis to let it lie on the z-zq plane pO| |, cf. 
Fig.l. It can be seen that in the rotated frame. Fig's. 4, 
5, 6 (c), the Fg{(p) increases much faster as the dimin- 
ishing of phase space scale than that in the thrust- minor 
frame, Fig's. 4, 5, 6 (&). As discussed in Ref. this 
is because the thrust-minor axis is basically determined 
by the first hard giuon emission and taking this axis as 
X axis to count the azimuthal angle ip will highly reduce 
the fluctuation of the direction of first hard gluon emis- 
sion. After rotation, the correlation between x axis and 
the direction of first hard gluon emission is relaxed and 
the full dynamical fluctuations are exhibited. 

We could also rotate the x axis around z for a random 
angle in each event The resulting Fq{ip) turn out 

to be the same as those in the rotated frame with x on 
the z-zq plane shown in Fig's. 4, 5, 6 (c). This confirms 
the cylindrical symmetry of the fluctuation in ip after the 
correlation with the thrust-minor is relaxed. 



VII. CONCLUSIONS 

The following conclusions could be drawn from the 
above investigation: 

1) The horizontal factorial moments (HFM) are equiv- 
alent to the vertical ones (VFM) only after the cumulant- 
variable transformation. Therefore, in the study of non- 
linear phenomena (intermittency or fractal) in high en- 
ergy collisions the HFM could be used only in combina- 
tion with the cumulant variables. On the other hand, the 
HFM is in its own right useful in single-event anaylsis. It 
can be seen from Eq.(4) that HFM is the average of the 
so called single-event factorial moments /g'^'' 

Hm) = (/i"V(/fV- (7) 

M 
m—l 



where rim is the number of particles of a single event 
in the mth sub-cell. The fluctuation of the single-event 
factorial moments around its average (HFM) is a charac- 
teristic of single-event fluctuations psl . 

2) The scaling properties of factorial moments in trans- 
verse directions {ip,pt) are very sensitive to the correction 
method used. They are unstable when using the correc- 
tion factor method, Eq.(5), even after the transformation 
Pt Inpt has been made. Using this method, the VFM 
in ip and the HFM in pt fall down instead of increase with 
the diminishing of phase space scale, while at the same 
time the HFM in ip and the VFM in pt do increase with 
the diminishing of scale, cf. Fig's. 4 (b) and (d). 

3) In the full phase space, the first few points of fac- 
torial moments do not lie on the scaling curve with the 
other points, due to the momentum conservation effect. 
This effect can be eliminated either through neglecting 
these points or through a cut in phase space. 

4) The thrust (or sphericity) major-minor frame is in- 
appropriate for the study of the scaling property of the 
azimuthal angle ip, because this frame is strongly cor- 
related with the direction of first hard giuon emission. 
Rotate the x axis to let it lie on the z-zq plane or rotate 
it randomly for each event can relax this correlation and 
exhibit the full dynamical fluctuations in ip. 

Therefore, the cumulant variables together with a 
frame rotated around the thrust (or sphericity) axis 
is the best for the investigation of the nonlinear phe- 
nomena (anomalous scaling of probability moments) in 
e+e^ collisions. The VFM and HFM are equivalent in 
this case. 
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